A Massive Data Parallel Computational Framework for Petascale/Exascale Hybrid Computer Systems

Heterogeneous systems are becoming more common on High Performance Computing (HPC) systems. Even using tools like CUDA and OpenCL it is a non-trivial task to obtain optimal performance on the GPU. Approaches to simplifying this task include Merge (a library based framework for heterogeneous multi-core systems), Zippy (a framework for parallel execution of codes on multiple GPUs), BSGP (a new programming language for general purpose computation on the GPU) and CUDA-lite (an enhancement to CUDA that transforms code based on annotations). In addition, efforts are underway to improve compiler tools for automatic parallelization and optimization of affine loop nests for GPUs and for automatic translation of OpenMP parallelized codes to CUDA. In this paper we present an alternative approach: a new computational framework for the development of massively data parallel scientific codes applications suitable for use on such petascale/exascale hybrid systems built upon the highly scalable Cactus framework. As the first non-trivial demonstration of its usefulness, we successfully developed a new 3D CFD code that achieves improved performance.Comment: Parallel Computing 2011 (ParCo2011), 30 August -- 2 September 2011, Ghent, Belgiu

From Physics Model to Results: An Optimizing Framework for Cross-Architecture Code Generation

Starting from a high-level problem description in terms of partial differential equations using abstract tensor notation, the Chemora framework discretizes, optimizes, and generates complete high performance codes for a wide range of compute architectures. Chemora extends the capabilities of Cactus, facilitating the usage of large-scale CPU/GPU systems in an efficient manner for complex applications, without low-level code tuning. Chemora achieves parallelism through MPI and multi-threading, combining OpenMP and CUDA. Optimizations include high-level code transformations, efficient loop traversal strategies, dynamically selected data and instruction cache usage strategies, and JIT compilation of GPU code tailored to the problem characteristics. The discretization is based on higher-order finite differences on multi-block domains. Chemora's capabilities are demonstrated by simulations of black hole collisions. This problem provides an acid test of the framework, as the Einstein equations contain hundreds of variables and thousands of terms.Comment: 18 pages, 4 figures, accepted for publication in Scientific Programmin

GeauxDock: A novel approach for mixed-resolution ligand docking using a descriptor-based force field

© 2015 Wiley Periodicals, Inc. Molecular docking is an important component of computer-aided drug discovery. In this communication, we describe GeauxDock, a new docking approach that builds on the ideas of ligand homology modeling. GeauxDock features a descriptor-based scoring function integrating evolutionary constraints with physics-based energy terms, a mixed-resolution molecular representation of protein-ligand complexes, and an efficient Monte Carlo sampling protocol. To drive docking simulations toward experimental conformations, the scoring function was carefully optimized to produce a correlation between the total pseudoenergy and the native-likeness of binding poses. Indeed, benchmarking calculations demonstrate that GeauxDock has a strong capacity to identify near-native conformations across docking trajectories with the area under receiver operating characteristics of 0.85. By excluding closely related templates, we show that GeauxDock maintains its accuracy at lower levels of homology through the increased contribution from physics-based energy terms compensating for weak evolutionary constraints. GeauxDock is available at http://www.institute.Loni.org/lasigma/package/dock/

A Massive Data Parallel Computational Framework for Petascale/Exascale Hybrid Computer Systems

large scale scientific applicatio

ITGB5 and AGFG1 variants are associated with severity of airway responsiveness

Background: Airway hyperresponsiveness (AHR), a primary characteristic of asthma, involves increased airway smooth muscle contractility in response to certain exposures. We sought to determine whether common genetic variants were associated with AHR severity. Methods: A genome-wide association study (GWAS) of AHR, quantified as the natural log of the dosage of methacholine causing a 20% drop in FEV1, was performed with 994 non-Hispanic white asthmatic subjects from three drug clinical trials: CAMP, CARE, and ACRN. Genotyping was performed on Affymetrix 6.0 arrays, and imputed data based on HapMap Phase 2, was used to measure the association of SNPs with AHR using a linear regression model. Replication of primary findings was attempted in 650 white subjects from DAG, and 3,354 white subjects from LHS. Evidence that the top SNPs were eQTL of their respective genes was sought using expression data available for 419 white CAMP subjects. Results: The top primary GWAS associations were in rs848788 (P-value 7.2E-07) and rs6731443 (P-value 2.5E-06), located within the ITGB5 and AGFG1 genes, respectively. The AGFG1 result replicated at a nominally significant level in one independent population (LHS P-value 0.012), and the SNP had a nominally significant unadjusted P-value (0.0067) for being an eQTL of AGFG1. Conclusions: Based on current knowledge of ITGB5 and AGFG1, our results suggest that variants within these genes may be involved in modulating AHR. Future functional studies are required to confirm that our associations represent true biologically significant findings

Lung eQTLs to Help Reveal the Molecular Underpinnings of Asthma

Genome-wide association studies (GWAS) have identified loci reproducibly associated with pulmonary diseases; however, the molecular mechanism underlying these associations are largely unknown. The objectives of this study were to discover genetic variants affecting gene expression in human lung tissue, to refine susceptibility loci for asthma identified in GWAS studies, and to use the genetics of gene expression and network analyses to find key molecular drivers of asthma. We performed a genome-wide search for expression quantitative trait loci (eQTL) in 1,111 human lung samples. The lung eQTL dataset was then used to inform asthma genetic studies reported in the literature. The top ranked lung eQTLs were integrated with the GWAS on asthma reported by the GABRIEL consortium to generate a Bayesian gene expression network for discovery of novel molecular pathways underpinning asthma. We detected 17,178 cis- and 593 trans- lung eQTLs, which can be used to explore the functional consequences of loci associated with lung diseases and traits. Some strong eQTLs are also asthma susceptibility loci. For example, rs3859192 on chr17q21 is robustly associated with the mRNA levels of GSDMA (P = 3.55 × 10(-151)). The genetic-gene expression network identified the SOCS3 pathway as one of the key drivers of asthma. The eQTLs and gene networks identified in this study are powerful tools for elucidating the causal mechanisms underlying pulmonary disease. This data resource offers much-needed support to pinpoint the causal genes and characterize the molecular function of gene variants associated with lung diseases

Survey of Period Variations of Superhumps in SU UMa-Type Dwarf Novae

We systematically surveyed period variations of superhumps in SU UMa-type dwarf novae based on newly obtained data and past publications. In many systems, the evolution of superhump period are found to be composed of three distinct stages: early evolutionary stage with a longer superhump period, middle stage with systematically varying periods, final stage with a shorter, stable superhump period. During the middle stage, many systems with superhump periods less than 0.08 d show positive period derivatives. Contrary to the earlier claim, we found no clear evidence for variation of period derivatives between superoutburst of the same object. We present an interpretation that the lengthening of the superhump period is a result of outward propagation of the eccentricity wave and is limited by the radius near the tidal truncation. We interpret that late stage superhumps are rejuvenized excitation of 3:1 resonance when the superhumps in the outer disk is effectively quenched. Many of WZ Sge-type dwarf novae showed long-enduring superhumps during the post-superoutburst stage having periods longer than those during the main superoutburst. The period derivatives in WZ Sge-type dwarf novae are found to be strongly correlated with the fractional superhump excess, or consequently, mass ratio. WZ Sge-type dwarf novae with a long-lasting rebrightening or with multiple rebrightenings tend to have smaller period derivatives and are excellent candidate for the systems around or after the period minimum of evolution of cataclysmic variables (abridged).Comment: 239 pages, 225 figures, PASJ accepte

Adaptation, spread and transmission of SARS-CoV-2 in farmed minks and associated humans in the Netherlands

In the first wave of the COVID-19 pandemic (April 2020), SARS-CoV-2 was detected in farmed minks and genomic sequencing was performed on mink farms and farm personnel. Here, we describe the outbreak and use sequence data with Bayesian phylodynamic methods to explore SARS-CoV-2 transmission in minks and humans on farms. High number of farm infections (68/126) in minks and farm workers (>50% of farms) were detected, with limited community spread. Three of five initial introductions of SARS-CoV-2 led to subsequent spread between mink farms until November 2020. Viruses belonging to the largest cluster acquired an amino acid substitution in the receptor binding domain of the Spike protein (position 486), evolved faster and spread longer and more widely. Movement of people and distance between farms were statistically significant predictors of virus dispersal between farms. Our study provides novel insights into SARS-CoV-2 transmission between mink farms and highlights the importance of combining genetic information with epidemiological information when investigating outbreaks at the animal-human interface

Shared genetic variants suggest common pathways in allergy and autoimmune diseases.

BACKGROUND: The relationship between allergy and autoimmune disorders is complex and poorly understood. OBJECTIVE: To investigate commonalities in genetic loci and pathways between allergy and autoimmune diseases to elucidate shared disease mechanisms. METHODS: We meta-analyzed two GWAS on self-reported allergy and sensitization comprising a total of 62,330 individuals. These results were used to calculate enrichment for SNPs previously associated with autoimmune diseases. Furthermore, we probed for enrichment within genetic pathways and of transcription factor binding sites, and characterized commonalities in the variant burden on tissue-specific regulatory sites by calculating the enrichment of allergy SNPs falling in gene regulatory regions in various cells using Encode Roadmap DHS data, and compared the allergy data with all known diseases. RESULTS: Among 290 loci previously associated with 16 autoimmune diseases, we found a significant enrichment of loci also associated with allergy (p=1.4e-17) encompassing 29 loci at a false discovery rate<0.05. Such enrichment seemed to be a general characteristic for all autoimmune diseases. Among the common loci, 48% had the same direction of effect for allergy and autoimmune diseases. Additionally, we observed an enrichment of allergy SNPs falling within immune pathways and regions of chromatin accessible in immune cells that was also represented in autoimmune diseases, but not in other diseases. CONCLUSION: We identified shared susceptibility loci and commonalities in pathways between allergy and autoimmune diseases, suggesting shared diseases mechanisms. Further studies of these shared genetic mechanisms might help understanding the complex relationship between these diseases, including the parallel increase in disease prevalence